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ABSTRACT 

We re-examine the stellar kinematics of the Solar neighbourhood in terms of the velocity 
Vq of the Sun with respect to the local standard of rest. We show that the classical determi- 
nation of its component V e in the direction of Galactic rotation via Stromberg's relation is 
undermined by the metallicity gradient in the disc, which introduces a correlation between 
the colour of a group of stars and the radial gradients of its properties. Comparing the local 
stellar kinematics to a chemodynamical model which accounts for these effects, we obtain 
(U, V, W) @ = (11.1+° ®, 12.24+°-^, 7.25+°-^) km s _1 , with additional systematic uncertainties 
~ (l,2,0.5)kms _1 . In particular, V Q is 7kms~ 1 larger than previously estimated. The new 
values of {U, V, W) are extremely insensitive to the metallicity gradient within the disc. 

Key words: stars: kinematics - Solar neighbourhood - Galaxy: fundamental parameters - 
Galaxy: kinematics and dynamics - Galaxy: disc 



1 INTRODUCTION 

The Sun's velocity v Q with respect to the Local Standard of Rest 
(LSRfl is required to transform any observed heliocentric velocity 
to a local galactic frame. Since this transformation is often neces- 
sary for scientific interpretation of observed velocities in terms of 
Galactic structure, the determination of v Q is a fundamental task 
of Galactic astronomy. The radial and vertical components U Q and 
Hq of v Q are straightforwardly obtained from the mean heliocentric 
velocities of several different groups of Solar-neighbourhood stars: 
t/ Q and Wo are simply the negative radial and vertical components 
of these mean^j. 

The component V Q of v Q in the direction of Galactic rotation 
is much harder to determine, because the mean lag with respect to 
the LSR, the asymmetric drift v a , depends on the velocity disper- 
sion a of the respective stellar population. The classical solution 
to this problem exploits the empirical linear relation between the 
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1 The LSR is the rest frame at the location of the Sun of a star that would 
be on a circular orbit in the gravitational potential one would obtain by az- 
imuthally averaging away non-axisymmetric features in the actual Galactic 
potential. 

2 According to the above definition, the LSR's radial and vertical motion 
w.r.t. the Galactic centre vanish. Therefore, the determination of f/ and 
W e from such means implicitly assumes that the Solar neighbourhood as 
a whole does not move radially or vertically w.r.t. the Galaxy. That such 
motions are at most sma ll is suggested by the proper motion of SgrA* 
teeid & Brunthaler 2004) and the mean radial velocity of the stars orbit- 
ing it (e.g. [geid et al. 2007). Moreover, such motions should also obey an 
asymmetric-drift like relation (see below), i.e. the mean velocities depend 
systematically on velocity dispersion, which is not observed. 



negative mean heliocent ric azimuthal vel ocity of any stellar sample 
v s = v- d + V Q and its rr 2 dStrom berg 1946). Hence, a straight-line fit 
yields V Q as the value of v s for a hypothetical population of stars on 
circular orbits, for which cr = 0. 

The theoreti cal underpinning of this me thod is the asymmetric 
drift relation (see lBinnev & Trema ine 2008, eq. 4.228) 
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where R is Galactocentric cylindrical radius, z the height above the 
plane, v c the circular speed, and v the number density of stars, while 
a bar indicates a v-weighted local mean. The equation applies sep- 
arately to each relaxed stellar population, for example to M stars 
or G stars. The idea behind the classical determination of V Q is that 
the square bracket in equation ([TJ takes essentially identical values 
for each stellar population, with the consequence that a plot of v s 
agai nst vj should be linear. 

iDehnen & Binnevl d 19981 . hereafter DB98) applied this method 
to a sample of ~ 15 000 main-sequence stars from the Hippar- 
cos catalogue and their value of V G = (5.25 ± 0.62) km s~' has 
been widely used. Recent re-det erminations us i ng an improved re- 
duction of the Hipparcos data ( Ivan Leeuwenl |2007|) confirm the 
DB98 value though wi th reduced error bars ( Ivan Leeuwenl 120071 : 
lAumer & Binnevlr2009l) . 

Ho wever, two rece nt studies call the DB98 value for V e into 
question. iBinnevI d2009l hereafter B09) fitted d istribution-function 
models (a) to velocity distributions inferred bv llvezic et alj ( 2008) 
from proper motions and photometric distances of stars in the Sloan 
Digital Sky Survey, and (b) to the spac e velocities of F and G 
in the Geneva-Copenhagen Survey (GCS. lNordstrom et alj|2004h . 
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The GCS stars are a subset of the Hipparcos stars (analysed by 
DB98) for which radial velocities have been obtained. B09 was 
able to obtain satisfactory fits to these data only if V a was larger 
than the DB98 value by ~ 6kms~', about ten times the formal 
error on V . Another body of evidence against the DB98 value 
for V B originates from radio-freq uency astrometry of masers in re- 
gions of massive-star formation (Rvgl et al. 2009;; lReid M. J. et al.l 
2009). If the DB98 value for V e is correct, th ese sources sys- 
temat ically lag circular rotation by ~ 17kms~' dReid M. J. et al.l 
|2009|). Such a high systematic lag is unexpected for young stars 
and lMcMillan & Binney! [2009) argued that a more plausible inter- 
pretation of the data is obtained if V Q exceeds the DB98 value by 
~ 6kms~'. 

This paper does two things: (i) it explains why the approach 
to the determination of V a by DB98 and subsequent studies is 
misleading, and (ii) it determines V Q from similar data but a dif- 
ferent methodology. Both these tasks are accomplished with the 
help of a particular chemo-d ynamical model of the Galaxy, that of 
ISchonrich & Bin ney (2009a, hereafter SB09a), but the points that 
we make are general ones and the role played by the SB09 model is 
essentially illustrative. In Section|2]we show that a metallicity gra- 
dient in the disc gives rise to distributions of mean azimuthal veloc- 
ity and velocity dispersion within the colour-magnitude plane that 
are much more complex than one naively expects, and we show that 
these distributions invalidate the methodology of DB98. In Sec- 
tion [3] we re-estimate V Q by fitting the entire velocity distribution 
of the GCS stars to the distribution predicted by the SB09 model 
without reference to the Stromberg relation. 



2 KINEMATICS IN COLOUR AND MAGNITUDE 

DB98 divided their sample of Hipparcos main-sequence stars into 
populations with different velocity dispersions by binning in B - V 
colour because colour is correlated with age and therefore with ve- 
locity dispersion. To examine the relation between colour, mean 
rotation velocity and velocity dispersion for stars near the Sun, 
we employ the SB09a model of the Galactic Disc. This model de- 
scribes the chemodynamical evolution of the thin a nd thick Galac- 
tic dis cs an d is a refinement of models pioneered bv lvan den Bergbl 
d 19621) and ISchmidtl dl963l) . The disc is divided into 80 annuli, 
within each of which the chemical composition of the ISM evolves 
in response to the ejection of material by dying stars, while stars 
form continuously with the current composition of the ISM. The 
new features of the model are (a) stochastic stellar accelerations ac- 
counting for heating processes; (b) radial stellar migration account- 
ing for both non-circular orbits and guiding-ce ntre shifts caused by 
stoch astic resonant scattering off spiral arms JSellwood & Binnevl 
2002); and (c) transfer of gas between annuli, both as result of res- 
onant scattering by spiral arms and as a result of a secular tendency 
of gas to spiral inwards through the disc. Surprisingly, the model 
contains both thin and thick di scs that are consistent with th e avail- 
able observational constraints dSchonrich & Binn ey 2009b). 

Fig-E shows the model kinematics in the colour-magnitude 
diagram. Each point in colour-magnitude space defines a separate 
sub-population whose asymmetric drift v a = v c - and radial ve- 
locity dispersion are plotted via colour coding, such that dynam- 
ically cold and warm populations are shown with blue and red 
shades, respectively. The region in the colour-magnitude diagram 
shown in this figure corresponds to the cuts used by DB98 to de- 
fine their sample. 

Our naive expectation is that as we proceed down the main 
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Figure 1. The variation with colour and magnitude of the asymmetric drift 
Ua (top) and radial velocity dispersion (bottom) in the SB09a model of the 
Solar neighbourhood. Shown is the range of colours and magnitudes used 
by DB98 to generate their main-sequence sample. Note that the number 
density of stars is highly non-uniform across the region shown. 



sequence from its blue end towards the main-sequence turnoff at 
B-V ~ 0.6, we encounter successively older stars with lower mean 
rotation velocities and higher velocity dispersions, so in both panels 
of Fig.Q]the shading should become redder as we move from left to 
right along the main sequence. The pattern actually found in Fig.Q] 
is more complex. Most notably, there is a pronounced velocity gra- 
dient across the lower main sequence. In the range 0.6 < B - V < 1 
the lower edge of the main sequence is dynamically warm (orange 
in Fig.QJ on account of subdwarfs, which are metal-poor and there- 
fore old with large velocity dispersions and low mean rotation rates. 
The number of these subdwarfs is small, however, so they will not 
have a significant impact on a sample binned by colour alone. More 
significant is the orange shading on the upper edge of the lower 
main sequence, which reflects the metallicity gradient within the 
disc: as metallicity increases, the main sequence shifts to the right, 
so in the upper panel the orange upper edge of the main sequence 
implies that the more metal-rich stars of the Solar neighbourhood 
are rotating more slowly because they formed at R < Ro. To the 
left of B-V - 0.5 this trend is weakened by contributions from 
old, sometimes metal-poor populations whose isochrones move up 
through this region. Still the more metal-rich main-sequence stars 



Local Kinematics and the Local Standard of Rest 3 





16 f 




14 




12 




10 






E 






8 


P 


6 




4 




2 





+ 


- 










+ 

+ + + 


+ + + + ' 




+ 


+ 

+ 




+ 


+ 


+ 

++++++++++++ 


- 






Figure 2. The asymmetric drift (top) and radial velocity dispersion (bot- 
tom) for stellar samples of given B-V colour drawn from the model thus 
simulating the effects of the Hipparcos and DB98 selection criteria. 



with smaller guiding centre radii give rise to slightly higher disper- 
sions and asymmetric drifts to the red side of the main sequence. 

The upper panel of Fig.[T]shows that in the crucial colour range 
0.4 < B-V < 0.6, the asymmetric drift is a complex function 
of colour and absolute magnitude because in this region stars of 
widely differing ages and metallicities are found as a result of old, 
metal-poor isochrones intersecting younger, metal-rich isochrones. 

The horizontal branch is clearly visible on both panels of 
Fig.Q]as an almost horizontal feature just below M v = 0. It is less 
pronounced in the upper panel because the blue end of the hori- 
zontal branch contains metal-poor stars, which tend to have large 
guiding-centre radii and therefore low v a even at large Vg. 

Fig.[2] shows the asymmetric drift u a and velocity dispersion 
obtained when stars are binned by colour alone. The lower panel 
can be compared with corre sponding observational plots in DB98 
and lAumer & Binnevl (2009). The model reproduces the structure 
of the data very well - in particular, the steepening in the slope 
around B-V = 0.4 and the flatness redwards of B-V = 0.6. 
The peak in velocity dispersion seen in the lower panel of Fig.|2] 
is much less evident in Fig. 2 of lAumer & Binnevl (2009) but can 
be traced in their cr R and <r- data. Note that the rise with B-V in 
for B-V < 0.4 is not accompanied by any change in v a . This 
unexpected phenomenon arises because at these colours the contri- 
bution of old metal-poor stars is increasing with B-V, and because 
we see many of these stars near pericentre, they have small asym- 
metric drifts despite their large random velocities. 

Since at its bright end the Hipparcos sample is close to being 
volume limited, the relative number of horizontal-branch stars is 
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Figure 3. Green squares: the asymmetric drift for synthetic stellar sub- 
samples defined by B-V colour plotted against their radial velocity disper- 
sion squared for the SB09a model. The red dots: the same relation obtained 
from a model with only one chemical composition (solar). Blue data points: 
the values of Aumer & Binney (2009) shifted by 11 kms~' and with the 
radial velocity dispersion increased by 7%. Purple line: a linear fit to the 
Hipparcos points in the range of t| used by DB98. 



small and the complex structure above the main sequence in Fig.Q] 
has no significant effect in Fig.|2] 

In Fig.[3] the squares show the resulting plot of the asymmetric 
drift u a against velocity dispersion squared for the synthetic sam- 
ples of Fig.|2j they do not lie on a straight line. The red dots show 
the plot one obtains if all stars are assigned solar metallicity. These 
dots do lie on a good approximation to a straight line0 The effect of 
re-assigning stars with large guiding-centre radii from solar metal- 
licity to their true, low metallicities is to move them from redder to 
bluer bins. Since these stars have small or even negative values of 
t/ a on account of the metallicity gradient in the disc, the transfer re- 
duces u d for the young, bluer bins and increases it in the old, redder 
bins. Consequently, the transfer morphs the near-straight line of the 
red points into the curve defined by the green squares. 

DB98 estimated V Q by fitting a straight line to the obser- 
vational analogue of Fig.[3] which is a plot of the solar veloc- 
ity relative to a colour-selected group of stars versus the squared 
velocity dispersion of that group. The blue data points in Fig.[3] 
show such data f or the Hipparcos sample in the re-analysis of 
lAumer & Binnevl l2009h after subtracting 1 1 km s 1 from each 
value of solar velocity. We see that for > 600 [km s~'] 2 the Hip- 
parcos data define the same straight line as the green crosses from 

the model. This straight line intercepts the v. d axis at 7kms~' 

rather than 0, causing V Q to be underestimated by this amount. For 
i% < 400 [kms -1 ] 2 the Hipparcos data points in Fig.[3]deviate from 
this straight line, but DB98 ignored samples with very low veloc- 
ity dispersion on the grounds that such samples may not be in dy- 
namical equilibrium. Indeed, both dissolving star clusters and the 
non-axisymmetric gravitational potentials of spiral arms are liable 
to distort the kinematics of stellar samples with low random veloc- 
ities such that equation (T) does not hold. 

The failure of the synthetic samples to follow a straight line 
in Fig.[3]implies that the square bracket in the asymmetric drift re- 



3 The slight deviation from a straight line of the model without metallicities 
is an effect of the approximations in SB09a and leads to a small underesti- 
mation of the real metallicity bias by the model. 
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Figure 4. Asymmetric drift velocity divided by squared velocity dispersion 
for synthetic colour-selected samples using the SB09a model. This ratio 
should be proportional to the square bracket in equation (T}. 



lation QJ does depend on colour: it varies by a factor ~ 4 in the 
colour range 0.4 < B - V < 0.6 as demonstrated in Fig.[4] A signif- 
icant contribution to the value of the bracket comes from the first 
derivative term, which is smallest for metal-poor populations be- 
cause their densities v decline more slowly outwards on account 
of the metallicity gradient in the disc. Physically, including metal- 
poor, thin-disc stars decreases u a because such stars typically visit 
the Solar neighbourhood at pericentre, where they have v$ > u c . 



3 DETERMINING THE SOLAR MOTION FROM THE 
VELOCITY DISTRIBUTION 

In view of the argument just presented, that the classical procedure 
cannot yield a reliable value of v Q , we now estimate v Q by fitting 
the observed distributions of heliocentric velocities to the veloc- 
ity distribution of the SB09a model. That is, we seek the offset 
—v Q from the circular velocity at which the model velocity dis- 
tributions provide the best match to the distribution of observed 
heliocentric velocities. B09 used an analogous procedure to argue 
that V^> a; 1 1 km s -1 ; however, the model distributions he used were 
obtained from analytic distribution functions rather from a model 
of the Galaxy's chemical evolution. By using velocity distributions 
that reflect much prior information about the chemodynamical his- 
tory of the Galaxy in place of simple analytic functions, we hope to 
achieve a closer fit to the observed velocity distributions and there- 
fore determine the requisite offset —v Q with greater precision. 

In Fig.|5]the points with ( Poisson) error bars are for a sub- 
sampl e of GC S stars for which lHolmberg. Nordstrom & Andersenl 
d2007l l2009h give reliable metallicities; thirty likely halo stars 
have been removed by requiring [Fe/H] > —1.2. This criterion is 
slightly stricter than that used in SB09b for the determination of 
the in-plane dispersion parameter (<r R of a 10 Gyr old local popula- 
tion), such that we now use a marginally smaller value, 43kms -1 , 
while lowering the vertical dispersion parameter to 23 km s -1 . The 
curves in Fig.[5]show the model distributions when offset by -v Q = 
(11.10+°;®, 12.24^, 7.25^) km s _1 - we used cubic splines to 
interpolate between individual data points provided by the model 
in the V component and determined the offset by maximising the 
likelihood of the data given the model. We used only five parame- 
ters for all three distributions: the two dispersion parameters and the 
three components of Solar motion, yet the fit is of good quality. The 
small fluctuations of the data around the model V and to a lesser 
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Figure 5. Curves: The model distributions predicted by the SB09a model in 
the U, V and W components of velocity (from top to bottom). Data points 
with Poisson error bars: the observed distributions of the GCS stars shifted 
by our estimate of Vq to optimise the fit of the data. 



extent U distributions are readily acc ounted for by the well known 
stellar moving groups toehneijl998h . likely caused by the dynami- 
cal influence of the Ga l actic bar and spiral structure jDehnenll999l : 
|Pe Simone et al j2004l : lAntoia et al.ll2009l : lMinchev et ai]|2009l) . 

A significant advantage of determining v Q from the entire 
sample, as in this section, rather than from subsamples as done 
in the past, is the robustness of the result to changes in the mod- 
elled metallicity gradient. In fact, eliminating the model's metallic- 
ity gradient changes V Q by less than 0. 1 km s -1 . 

One should note that the quoted errors on the components of 
v Q are purely formal. Sources of additional systematic error include 
the possible presence of halo stars in the sample, the dynamical 
approximations used in constructing the model, and the effects of 
stellar streams, which have a big impact on the observed distribu- 
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tion of stars near the circular velocity but are completely excluded 
from the model. Fortunately the likelihood of the data used here is 
not particularly sensitive to the fit of the data to the model around 
the peak density. In view of these uncertainties we roughly estimate 
systematic errors of ~ (1, 2, 0.5) km s _1 , assuming that U and W are 
mostly affected by distortions by streams and V showing even more 
structure and having an additional uncertainty from the modelling. 
This is in perfect agreement with previous estimates as regards W e 
and slightly higher in U Q compared to the DB98 value, which can 
be traced back of the larger influence of the Hercules stream at 

30kms~' (lowering the estimate) on their statistics. However, 

it differs significantly from the value for V n a* 5.2 ± 0.5 km s -1 ob - 
tained by the classical technique (DB98, lAumer & Binnevi r2009). 
Our value for is in good agreement with V - 11 kms~' pro- 
posed by B09. Given the residual uncertainties of v Q , it is question- 
able whether the standard practice of "correcting" observed (helio- 
centric) velocities for the Solar motion is useful, at least the adopted 
value should be explicitly provided. 



4 CONCLUSIONS 

The metallicity gradient in the Galactic disc causes a systematic 
shift in the kinematics especially near the turnoff region. By the re- 
lationship between the colour and metallicity of a star, the more 
metal-rich populations, with on average lower angular momen- 
tum and thus higher asymmetric drifts, are displaced relative to 
their metal-poor counterparts, which have lower asymmetric drifts. 
When stars are binned by colour, the metallicity gradient in the 
Galactic disc prevents the relationship between mean rotation speed 
and squared velocity dispersion taking the linear form predicted 
by a naive application of the Strbmberg relation. This breakdown 
in the conventionally assumed linearity invalidates the traditional 
technique for determining the Sun's velocity with respect to the 
LSR, which involves a linear extrapolation to zero velocity dis- 
persion of the empirical relation between the mean velocity and 
squared velocity dispersion. Moreover the SB09a model predicts 
that a treacherous linear relationship underestimating the solar az- 
imuthal motion by ~ 7kms~' will be mimicked redwards of the 
onset of the turnoff region, coinciding well with the behaviour ob- 
served in the Hipparcos data. 

The Sun's velocity with respect to the LSR may be 
alternatively determined from the velocity offset that opti- 
mises a model fit to the observed velocity distribution. Us- 
ing the velocity distribution predicted by the SB09a model of 
the chemodynamical evolution of the Galaxy, we find v Q = 
(ll.l!o?5. 12 - 24 -oat 7.25*0 36) km s" 1 and roughly estimate the 
systematic uncertainties as ~ (1, 2, 0.5) km s _I . The radial and verti- 
cal components of this value of v Q agree with earlier estimates, but 
the V component is larger than the widely used value of DB98 by 
~ 7 km s -1 . This is in nice concordance with the model expectations 
for the systematic error arising from naively using the Strbmberg 
relation and in good agreement with the value obtained by B09 us- 
ing a similar method but with a less sophisticated distribution func- 
tion. Curiou sly it agr e es we ll with the result V Q ~ 10 - 13kms~' 
obtained by IPelhaveN 19651) using the classical method with pre- 
Hipparcos data. 

In this paper we have relied heavily on the SB09a model, so 
the question arises of how vulnerable our argument is to the model's 
shortcomings. Our critique of the classical approach to the determi- 
nation of v Q is secure so long as the disc has significant age and/or 
metallicity gradients. It is beyond question that such gradients ex- 



ist, so the classical technique is certainly unreliable. Our proposed 
value of Vq is essentially independent of the assumed metallicity 
gradient, but does have some sensitivity to the dynamical approxi- 
mations used in making the SB09a model - plausible variations of 
how one handles the secular acceleration of stars lead to changes 
in the estimated value of V Q by 0.5 to 1 km s -1 . Modest reassurance 
that the error of our value of V a is less than 2kms~' is furnished 
by the fact that B09 favoured the same value using a distribution 
function that takes no account of the age and metallicity gradients 
in the disc. Consequently, our result is probably not sensitive to 
the assumptions about star formation and chemical evolution made 
by SB09a. However, as we develop more elaborate models of the 
Galaxy which fit a wider range of data, in particular more distant 
stars, we anticipate further small revisions in the value of V a . 
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